First analysis of anisotropic flow with Lee— Yang zeroes 
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We report on the first analysis of directed and elliptic flow with the new method of Lee- Yang 
zeroes. Experimental data are presented for Ru+Ru reactions at 1.69^4 GeV measured with the 
FOPI detector at SIS/GSI. The results obtained with several methods, based on the event- plane 
reconstruction, on Lee- Yang zeroes, and on multi-particle cumulants (up to 5th order) applied for the 
first time at SIS energies, are compared. They show conclusive evidence that azimuthal correlations 
between nucleons and composite particles at this energy are largely dominated by anisotropic flow. 

PACS numbers: 25.70.-z, 25.75.Ld 



The study of collective flow in relativistic heavy ion 
reactions is of great interest since it is expected to shed 
light on our knowledge about the properties of hot and 
dense nuclear matter and the underlying equation of state 
(EoS) 0. As pointed out early on, nuclear collective flow 
is also influenced by the momentum-dependent interac- 
tions and the in-medium nucleon-nucleon cross section 
0,0. Both effects play a crucial role in the determina- 
tion of the EoS and cannot be neglected at intermediate 
energies. In this regard both directed and elliptic flow are 
a field of intense experimental and theoretical researches 
(see Q and references therein). 

Most flow analyses, based cither on the reaction plane 
reconstruction (the so-called event-plane method) [5| or 
on two-particle azimuthal correlations [f| rely on the as- 
sumption that the only correlations are those stemming 
from the existence of the reaction plane. Other correla- 
tions (usually called non-flow) , such as small-angle corre- 
lations due to final state interactions and quantum sta- 
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tistical effects 0, correlations due to resonance decays 
y| and mini-jet production are neglected. In recent 
years, several alternative techniques were introduced, in 
which non-flow correlations can be unraveled. The cumu- 
lant method is based on a cumulant expansion of multi- 
particle (typically four particles) correlations fT(i| , which 
eliminates most of non-flow correlations. It has been ap- 
plied at ultra-relativistic energies, at RHIC and SPS for 
directed and elliptic flow studies and also for higher har- 
monic measurements 0,0,0,0,^3- More recently, a 
new method based on an analogy with the Lee- Yang the- 
ory of phase transitions |16(, where flow is extracted di- 
rectly from the genuine correlation between a larg e num- 
ber of particles, has been proposed 0, IlsL Il9| This 
method is expected to provide the cleanest separation 
between flow and non-flow effects. 

We present the first analysis of collective flow using the 
new method of Lee-Yang zeroes. The cumulant method 
is also applied, for the first time at SIS energies. A 
comparison with results obtained with the event-plane 
method is performed. We are thus able to check for the 
first time the validity of standard methods at SIS ener- 
gies, by investigating possible contributions of correla- 
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tions unrelated to the reaction plane which could intro- 
duce distortions on directed and elliptic flow results. 

The data set presented in this work concerns Ru+Ru 
reactions at 1.69^4 GeV measured with the FOPI de- 
tector installed at the SIS accelerator facility of GSI- 
Darmstadt. FOPI is an azimuthally symmetric appara- 
tus made of several sub-detectors which provide charge 
and mass determination over nearly the full 47r solid an- 
gle. The central part (33° < #i a b < 150°) is placed in a 
super-conducting solenoid and consists of a drift cham- 
ber (CDC) surrounded by a barrel of plastic scintilla- 
tors. Particles measured in the CDC are identified by 
their mass using magnetic rigidity and energy loss. The 
forward part is composed of a wall of plastic scintilla- 
tors (1.2° < #i a b < 30°) and an other drift chamber 
(Helitron) mounted inside the super-conducting solenoid. 
The plastic wall provides charge identification of the re- 
action products, combining time of flight and specific en- 
ergy loss informations. For the present analysis, the for- 
ward wall and the CDC were used. More details on the 
configuration and performances of the different compo- 
nents of the FOPI apparatus can be found in [2(j • 

The events are sorted out according to their degree 
of centrality by imposing conditions on the multiplicity 
of charged particles measured in the outer plastic wall 
(7° < 0i ab < 30°) H3, named PMUL. The flow analy- 
sis presented here was carried out for about 2.9 million 
events belonging to the centrality class labelled PMUL4, 
which corresponds to a mean geometrical impact param- 
eter of 2.9 fm and to a geometrical impact parameter 
range from 1.6 fm to 3.9 fm, obtained assuming a sharp- 
cut-off approximation [22| . 

We recall that directed flow («i) and elliptic flow {V2) 
are quantified by Fourier coefficients of the azimuthal dis- 
tributions [23^. v n — (cosn(tp — (fin)), where ip is the 
particle azimuthal angle and ip R is the azimuth of the 
reaction plane. 

In the conventional method, the reaction plane is es- 
timated event by event according to the standard trans- 
verse momentum procedure devised in p|, which allows 
to construct the event-plane vector 
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(i) 



The sum runs over all charged particles in the event, ex- 
cept pions identified in the CDC. u„ is the unit vector 
parallel to the particle transverse momentum (i.e. u„ = 
(cos (ft, , sin (/?„), where (/?„ is the particle azimuth), and u> v 
is a weight to improve the resolution, depending on the 
scaled center-of-mass (cm.) rapidity y(°) = (y/y p ) c .m. 
(the subscript p refers to the projectile): ld v = — 1 for 
< -0.3, uj v = +1 for yW > 0.3, and uj v = other- 
wise. The azimuth of Q, denoted by ^f R , is an estimate 

of CfR. 

The Fourier coefficients v n are calculated using the for- 
mula 



v n {EP} 



{cosn(ip - V R )) 
(cosnAip R ) 



where {EP} stands for "event-plane" . In the numerator, 
the particle of interest is excluded from the sum in Eq. 
to avoid autocorrelation effects. The resolution factor 
1/ (cosnAif R ) is an estimate of the error A(p R — fy R — (p R 
on the determination of the reaction plane. This factor is 
calculated according to the procedure proposed in [24J, 
and involves the correlation between randomly chosen 
sub-events [j| . The numerical values for the PMUL4 cen- 
trality class are 1/ (cos Aip R ) = 1.17 for directed flow and 
l/(cos2A^7?) = 1.69 for elliptic flow, corresponding to a 
resolution parameter \ ~ 1-47 |24j . Several procedures 
have been developed in order to take into account correla- 
tions due to overall transverse momentum conservation, 
by using standard methods [2^, |2(| . Here, to subtract 
these correlations the event-plane method has been im- 
proved by introducing a recoil correction, as proposed in 

Let us now recall the principle of the Lee- Yang zeroes 
procedure to analyze flow. A more complete description 
of the method can be found in 0, 0, 0] . The method 
is based on the location of the zeroes, in the complex 
plane, of a generating function of azimuthal correlations, 
in close analogy with the theory of phase transitions of 
Lee and Yang [Tg. The first step of the procedure is to 
determine the "integrated" directed flow, defined as the 
average projection of Q on the (true) reaction plane 



Vi = (Q x cos ip R + Q v sin ip R ) e 



(3) 



where Q x , Q y are the components of Q, and the average 
is taken over events in a centrality class. For this purpose, 
one introduces the following complex-valued generating 
function [HSf : 



G e {ir) = / J|[l + iru v cos(^ - 9)} 



(4) 



where r is a positive real variable, 9 is an arbitrary ref- 
erence angle, uj v is the same weight as in Eq. QJ, and tp v 
is the particle azimuthal angle. 

Figure^displays the amplitude of the generating func- 
tion |£r (w")| versus r for 9 = (circles). It starts at 
a value of 1 for r = and quickly decreases as r in- 
creases, which means that particles are strongly corre- 
lated: for uncorrelated particles, indeed, \G e (ir)\ is iden- 
tically equal to unity within statistical fluctuations. A 
sharp minimum of \G e (ir)\ then occurs, which is in fact 
compatible, within statistical fluctuations, with a zero of 
G (ir) (see insert in Fig.[TJ. Following the general argu- 
ments presented in [T7I Il8t IT ^ |. this is a clear indication 
that correlations are due to collective flow. The posi- 
tion Tq of the first minimum yields an estimate Vf of the 
integrated flow V\ : 



Tr8 _ 301 

v i ~ ~~F> 

7n 



(5) 



(2) 



where joi = 2.40483 is the first root of the Bessel func- 
tion Jo(x). Strictly speaking, like other flow analysis 
methods, the present one is only able to determine the 
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FIG. 1: \G e {ir)\ versus r for 9 = 0. A zoom of \G e (ir)\ 
(Eq. gj) around its first minimum is shown in the insert and 
the solid line is just to guide the eye. See text for details. 

absolute value of V\ . The sign is assumed to be positive 
at these energies. 

A potential limitation of the method comes from sta- 
tistical errors, which can be much larger than with the 
event-plane method. The reason why statistical errors 
depend on the method used is that the reaction plane 
is unknown, and that the v n are obtained through the 
indirect observation of a correlated emission p]. The 
statistical errors depend on the observables used to char- 
acterize this correlation. The important quantity here 
is the resolution parameter x, related to the well-known 
event-plane resolution [24[. If % > 1, which means that 
the reaction plane can be reconstructed with reasonable 
accuracy, all methods yield statistical errors of the same 
order of magnitude as if the reaction plane was exactly 
known while systematic errors from non-flow effects are 
expected to be much smaller with Lee Yang zeroes. If 
X < 0.5, on the other hand, statistical errors prevent the 
use of Lee- Yang zeroes. For the present analysis, we find 
X ~ 1.45, which definitely indicates that statistical errors 
are not a problem here (see Fig. [21 and Fig.^. 

An alternative form of the generating function, which 
can be used instead of Eq. is |l7j : 

G e (ir) — / exp ( ir u v co$((p v — 6) I ) . (6) 

\ V v / I events 

This gives the squares in Fig. ^ There is no obvious 
relation between Eq. (JIJ and Eq. ©. In particular, the 
latter differs from unity for r > even if particles are un- 
corrected because of "autocorrelation" terms. But quite 
remarkably, the first minimum occurs at the same place 
with either form. This is a further indication that it is 
due to flow, as anticipated in Ref. ( see m partic- 
ular Appendix A). The analysis of differential directed 
and elliptic flow presented below was carried out using 
Eq 0. The results obtained with Eq. JSJ) are very simi- 



lar, thereby confirming that the method is insensitive to 
autocorrelations. 

Once the first minimum Tq has been determined, the 
Fourier coefficients are estimated from the following 
equation: 

v e n <xRe/cosn(tp - 9)]J'[1 + ir^0J v cos(ip v - 6)]\ (7) 

where <p is the azimuthal angle of the analyzed particle, 
and the notation F| ' means that the particle of interest 
is excluded from the product in order to avoid autocor- 
relations. The average is over a particle type in a given 
phase-space region, in all events. A proportionality con- 
stant ensures that the result is consistent with the esti- 
mate of the i nteg rated flow Vf. Its expression can be 
found in Ref. . 

The procedure is repeated for several values of 9 (typ- 
ically, 5 equally spaced values from to 4-7r/5), and the 
results are found to be independent of 9 except for sta- 
tistical fluctuations. This demonstrates that the results 
are not affected by detector azimuthal asymmetries. The 
final estimates shown below are averaged over 9, which 
reduces the statistical errors by about a factor of 2. 

Before we come to the results, let us say a few words 
about the cumulant method 0,which has been already 
applied by several experiments [ll], [Tj, EiJ EH Ell ■ This 
method makes use of multi-particle correlations to esti- 
mate directed and elliptic flow. One can construct several 
independent estimates of v\ and i>2, depending on how 
many particles are correlated: 2 or 4 for v%, 3 or 5 for 
V2 ■ The number of particles involved is referred to as the 
order of the cumulant. Lowest-order estimates of v\ and 
V2 are not corrected for non-flow effects, and therefore 
are excepted to be similar to estimates from the event- 
plane method without recoil correction. The higher the 
order, the smaller the bias from non-flow correlations. 
Lee- Yang zeroes are essentially the limit of cumulants 
when the order goes to infinity, and therefore minimize 
the bias from non-flow effects. 

The features of differential directed and elliptic flow at 
SIS energies have been discussed in several publications 
[27L l2Sj. Here, we focus on the comparison between the 
different procedures investigated in this work. 

In the following figures only statistical errors are 
shown. Possible sources of systematic uncertainties have 
been studied using IQMD events [2^| passed through a 
complete GEANT simulation of the detector. We found 
that the full simulation underestimates v± of protons 
(deuterons) by about 6% (4%), in relative value, in the 
phase space region under consideration and for data inte- 
grated over transverse momentum (pt)- These distortions 
are mainly due to a track-density effect which leads to a 
loss of particles in the directed flow direction. They are 
independent of the procedure. 

The differential directed flow calculated with the 
method of Lee- Yang zeroes (circles) is shown in Fig. [5] 
for protons (upper panel) and deuterons (lower panel) in 
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FIG. 2: «i versus transverse momentum for protons (upper 
panel) and deuterons (lower panel) measured in semi-central 
events and in a rapidity window in backward hemisphere. See 
text for details. 



a rapidity window in the backward hemisphere. The val- 
ues are compared to those obtained from the event-plane 
analysis (squares). Also shown are the second (stars) and 
fourth (crosses) particle cumulant values. 

A first look at Fig. [21 shows that all methods give sim- 
ilar results. This proves that azimuthal correlations be- 
tween nucleons and composite particles at SIS energies 
are dominated by anisotropic flow and that non-flow cor- 
relations, if any, are of smaller magnitude. Figure Halso 
shows that cumulants and Lee- Yang zeroes can be suc- 
cessfully used to analyze anisotropic flow at SIS energies. 

A detailed examination of the results in Fig. |3 re- 
veals however that there are small differences between the 



methods, beyond statistical errors at high p t . First, there 
is a small difference between the event-plane method and 
the second-order cumulant. This is due to the recoil cor- 
rection for overall momentum conservation, which is ap- 
plied in the event-plane method, but not in the cumulant 
method. We have checked that the event-plane method 
without recoil correction and the second-order cumulant 
give compatible results. It is important to emphasize 
that four-particle cumulants and event-plane results dif- 
fer. The difference between second-order cumulant and 
fourth-order cumulant is also observed in analyses of el- 
liptic flow at RHIC ^3 ' where discrepancies are larger in 
relative value. There, it was suggested that the difference 
may be due to fluctuations of the flow within the sam- 
ple of events, corresponding to variations in the impact 
parameter or in the initial conditions [ill l3(l | . This ef- 
fect leads to smaller flow estimates with the four-particle 
cumulant than with the second-order cumulant, in abso- 
lute value, independently of p t . The opposite behavior is 
evidenced in Fig. |3 from which one concludes that such 
fluctuations are not responsible for the observed differ- 
ences. It seems therefore more likely that they are due to 
non-flow correlations. The fact that they increase with 
Pt suggests that they are mostly due to overall trans- 
verse momentum conservation |3l) . The recoil correction 
which has been used in the event-plane method to correct 
for this effect relies on a non-relativistic formalism |2^| 
and that may explain the difference relative to the fourth- 
order cumulant. Moreover, it is worth to stress the fact 
that overall momentum conservation, which is a long- 
range effect involving all par ticles, effectively behaves as 
a short-range correlation [171 l31j . As a consequence it is 
eliminated by using fourth-order cumulant. On the other 
hand, results from four-particle cumulants and Lee- Yang 
zeroes are perfectly compatible. This lends support to 
the idea that both methods are able to extract reliably 
the genuine collective flow at SIS energies. 

Figure [21 displays the proton differential elliptic flow, 
estimated from the same methods. The results now con- 
cern the mid-rapidity region (—0.2 < < 0.) and the 
corresponding p t range of the CDC acceptance. Effects 
of non-flow correlations such as correlations due to mo- 
mentum conservation are expected to be less pronounced 
on elliptic flow than on directed flow. Indeed, the differ- 
ences between the methods are smaller, in absolute value, 
than for directed flow and can be considered as almost 
negligible within statistical error bars. However the gen- 
eral trend seems to be that the lowest-order cumulant 
(stars) gives a slightly larger signal than Lee- Yang ze- 
roes (circles) and fifth-order cumulant (crosses) , at all pt ■ 
This small difference could be due to impact parameter 
fluctuations [Tj within the PMUL4 centrality bin which 
increase the estimates of vi from cumulants, in absolute 
value, independently of p t . This bias is expected to be 
more pronounced for the third-order cumulant than for 
high-order cumulants. v-i values from the event-plane 
method (squares) are not distorted by such fluctuations 
because of the high accuracy on the reaction plane de- 
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FIG. 3: «2 versus transverse momentum for protons measured 
in semi-central events and around mid-rapidity. See text for 
details. 



termination. Similar trends have been also obtained for 
deuterons. 

In summary, we have presented the first analysis of 



directed and elliptic flow in heavy-ion collisions using the 
method of Lee- Yang zeroes. Results were obtained from 
the FOPI experiment at GSI. Such method is expected to 
provide the best possible separation between correlations 
due to flow and other correlations. We were thus able to 
check explicitly that most azimuthal correlations between 
protons and composite particles at SIS energies are due to 
their correlation with the reaction plane of the collision. 
There is no evidence for event-by-event fluctuations of 
directed flow. Non-flow effects are small; they are clearly 
seen only on directed flow at high p t , and may be entirely 
ascribed to global transverse momentum conservation. 
They are eliminated using four-particle cumulants or the 
Lee- Yang zeroes procedure. Results were presented only 
for semi-central events, for sake of brevity. The analysis 
was also carried out for other centrality classes, covering 
an impact parameter range up to 7 fm, and led to similar 
conclusions. 

Such analysis is promising for studying pion flow. 
Since most pions originate from A decays, the resulting 
non-flow correlations with protons may contaminate 
the flow analysis if standard procedures are applied. 
The Lee- Yang zeroes method, which is insensitive to 
non-flow effects, should be able to provide reliable results. 
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